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Monte Carlo simulations are frequently utilized in radiation dose assessments. However, many researchers 
find the prevailing computing platforms to be intricate. This highlights a pressing need for a specialized frame- 
work for phantom dose evaluation. To address this gap, we developed a user-friendly radiation dose assessment 
platform using the Monte Carlo toolkit, Geant4. The Tsinghua University Phantom Dose (THUDosepp) aug- 
ments the flexibility of Monte Carlo simulations in dosimetric research. Originating from THUDose, a code 
with generic, functional, and application layers, THUDosepp focuses predominantly on anatomical phantom 
dose assessment. Additionally, it enables medical exposure simulation, intricate geometry creation, and supports 
both three-dimensional radiation dose analysis and phantom format transformations. The system operates on a 
multi-threaded parallel CPU architecture, with some modules enhanced for GPU parallel computing. Bench- 
mark tests on the ICRP reference male illustrated the capabilities of THUDosepp in phantom dose assessment, 
covering the effective dose, three-dimensional dose distribution, and three-dimensional organ dose. We also 
conducted a voxelization conversion on the polygon mesh phantom, demonstrating the method’s efficiency and 
consistency. Extended applications based on THUDosepp further underline its broad adaptability. This intu- 


itive, three-dimensional platform stands out as a valuable tool for phantom radiation dosimetry research. 


Keywords: Monte Carlo, phantom, dose assessment, voxelization, three dimensional 


I. INTRODUCTION 


With the growing adoption of nuclear techniques in ra- 
diation protection and medical physics[1—6], accurate dose 
calculations for human phantoms have become indispens- 
able. Phantom-dose evaluations are pivotal in clinical prac- 
tices such as radiotherapy, CT irradiation, X-ray radiography, 
and in radiation protection undertakings such as retrospec- 
tive dose reconstruction and epidemiological studies [7-9]. 
Among the myriad techniques available, the Monte Carlo al- 
gorithm is esteemed as the gold standard for dose compu- 
tation. This method encapsulates all physical processes and 
offers insights into transport processes by emulating particle- 
matter interactions[10]. Furthermore, the Monte Carlo ap- 
proach facilitates the translation of microscopic mechanisms 
into macroscopic doses and exhibits higher precision over an- 
alytical algorithms, especially in heterogeneous geometries. 

Many specialized Monte Carlo software tools are dedi- 
cated to phantom dose calculations, including GATE[11, 12], 
TOPAS[13, 14], and ARCHER[15-—17]. Each has been metic- 
ulously researched and optimized for specific application do- 
mains. Other notable methods in this arena are MCNPX[18], 
FLUKA[19], and Geant4[20]. While MCNPX and FLUKA 
largely operate via input card simulations, presenting a con- 
siderable workload for intricate geometries and source mod- 
eling, Geant4 stands out. Geant4 can accommodate more 
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complex geometries and offers an extensive suite of func- 
tionalities, enhancing code flexibility. Written in C++, this 
open-source platform facilitates easier packaging and further 
development. 

Geant4’s flexibility leads to complexity. Writing and de- 
bugging code for a novel simulation requires considerable ef- 
fort, and it can be difficult for researchers without significant 
programming experience. Additionally, although Geant4 sup- 
ports phantom dose calculations, the user should import the 
geometry, dose algorithm, and tally score. If Geant4 can be 
packaged, then a specific function can be connected to the in- 
terface parameters, and a convenient and stable phantom dose 
calculation platform can be generated for specific phantom 
dose calculations. 

Therefore, we developed a dedicated phantom dose calcu- 
lation platform, THUDosepp, based on the Geant4 toolkit, 
which includes the following: 


e preserve the underlying Geant4 code 
e visual user interface 
e provide various usage of phantom dose calculations 


e support CPU and GPU multi-threaded parallel comput- 
ing 


e support three-dimensional visualization of geometry 
and results 


The platform THUDosepp is tailored for pragmatic appli- 
cations in clinical medicine and radiation protection. It boasts 
a user-friendly interface that facilitates complex geometric 
modeling and supports the transport of electrons, photons, 
neutrons, protons, and heavy ions. Additionally, the platform 
accommodates both single and dual phantom dose calcula- 
tions, catering to the routine needs of medical physics and 
radiation protection. It offers methodologies for evaluating 
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Fig. 1. The overall structure of THUDose platform. 


organ doses from both external and internal exposures. Keep- 
ing in sync with the evolution of phantoms, THUDosepp 
has implemented a procedure for the voxelization of polygon- 
mesh phantoms. The platform further enhances user experi- 
ence by supporting three-dimensional visualization, spanning 
from the initial modeling phase to the final output. Grounded 
in a robust software architecture, THUDosepp has proven its 
merit in multiple research endeavors, covering scenarios of 
utmost interest to the majority of users and consistently offer- 
ing convenience and stability. 


I. MATERIAL AND METHODS 


The THUDosepp platform is based on THUDose[21], and 
was then expanded and developed for phantom dose calcula- 
tion and three-dimensional visualization. The extended pro- 
gram has additional features, such as creating and converting 
phantoms, pre-defining standard source types, external and 
internal irradiation dose algorithms, and geometry visualiza- 
tion. 


A. The universal MC program THUDose 


THUDose is an in-house general Monte Carlo procedure 
written in C++ that wraps and extends the Geant4 toolkit. 
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Fig. 2. The interactive interface organization structure of THUDose. 


The platform operates via a Graphical User Interface (GUI) 
through which projects are built. In the background, XML 
files were saved using the settings set by the interface. The 
general, function, and application layers are the three layers 
that constitute THUDose. Figure 1 illustrates the system’s 
overall structure. The four steps comprise a simulation in 
THUDose: geometric data, source definition, tally scores, 
and run parameters. The accompanying Figure 2 shows 
the organizational structure of the interactive interface and a 
schematic. 


In THUDose, the developers integrated physics lists de- 
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rived from the cross-sectional libraries of both Geant4 and 
ENDPF, selecting hadron and electromagnetic physics models 
for registration. This program was used to assess the trans- 
port of electrons, photons, protons, and neutrons. For material 
construction optimization, the platform provides two distinct 
strategies: The first approach utilizes the "Periodic Table" in- 
terface, allowing users to define materials on an individual 
basis. As an alternative, users can employ an Excel spread- 
sheet to craft a bill of materials, which facilitates the addition 
of materials in batches. 

The source definition in THUDose is versatile, accommo- 
dating any particle type, energy, orientation, or spatial dis- 
tribution. It allows the creation of multiple tallies for a sin- 
gle project, with each tally functioning independently. The 
software diligently logs energy deposition, flux, and particle 
count in the medium, segmented by particle type. To boost 
computational efficiency, the developers integrated parallel 
computing modules based on MPI. At present, THUDose can 
utilize up to 2000 parallel computing cores, contingent upon 
the computer’s memory capacity. Furthermore, the software 
is equipped with GPU parallel computing functionalities for 
electrons and photons[8, 22]. The team plans to introduce 
support for GPU-based proton transport in the near future. 


B. Extensions for phantom dose assessment 


Four processes comprise the majority of MC-based 
phantom dose assessments: phantom construction, parti- 
cle transportation, dose calculation algorithm, and visu- 
alization. We investigated these four areas and devel- 
oped a three-dimensional phantom dose assessment platform, 
THUDosepp based on THUDose. 


1. Phantom construction and conversion 


Currently, the voxel phantom stands as the predominant hu- 
man anatomy model employed in radiation dose assessment. 
For precise dose assessment in intricate workplaces, it is vital 
to produce phantoms that offer adjustable sizes and postures. 
Such malleable models can be crafted using the Boundary 
Representation (BREP) geometric definition methods, which 
can take the form of nonuniform rational B-splines (NURBS), 
polygon meshes, or a combination thereof. Relative to voxel 
models, the BREP model offers a broader range of deforma- 
tion operations. It can represent new anatomical structures for 
staff or patients even in the absence of comprehensive whole- 
body image data[23]. For the purpose of this work, all BREP 
models, which are expressed as polygon-mesh phantoms be- 
fore final computations, are termed ’ polygon-mesh phantom’. 

Direct calculations on polygon mesh phantoms are unfea- 
sible. Typically, they are converted into either voxels or 
tetrahedrons to facilitate particle transport[24, 25]. As of 
now, THUDosepp supports two such conversion techniques. 
However, given that tetrahedralization is not the focal point 
of this study, only the voxelization process is elucidated. 
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For the voxel phantom, THUDosepp can import data in a 
binary file format. Pertinent files include both a configuration 
file and binary file. Details, such as phantom’s size, voxel’s 
resolution, and phantom’s central coordinates, are specified 
in the configuration file. Each voxel’s value corresponds to a 
material index within a predefined list, indicating its material 
classification. 

For polygon-mesh phantoms, THUDosepp offers a trans- 
formation procedure, termed voxelization. This method con- 
verts a three-dimensional object’s geometric representation 
into a voxel format that closely mirrors the original. The 
resulting data captures both the surface and the inner geo- 
metric details. For rapid voxelization, THUDosepp employs 
OpenGL and utilizes the fragment shader in the rendering 
process to transition the intersection of triangles and voxels 
into sampling points. In parallel, computing techniques expe- 
dite the voxelization outcomes. 

To enhance clinical efficacy, precise radiation dose evalua- 
tions based on actual patient CT scans are paramount. To this 
end, a phantom construction module was developed to load 
DICOM files, extract HU values, and subsequently transform 
them into densities and materials for ensuing calculations. 


2. Particle transportation and source definition 


The THUDosepp physical process is based on the THU- 
Dose physical model, facilitating the simulation and transport 
of photons, electrons, neutrons, protons, and heavy ions. The 
program has predefined the eight irradiation types highlighted 
in ICRU Publication 95: antero-posterior (AP), posterior- 
anterior (PA), left lateral (LLAT), right lateral (RLAT), ro- 
tational (ROT), isotropic (ISO), superior hemisphere semi- 
isotropic (SS-ISO), and inferior hemisphere semi-isotropic 
(IS-ISO). For X-ray radiography, which is prevalent in med- 
ical applications, the source has been parameterized. Users 
simply define the source axis distance, deflection angle, and 
field size to initiate vertebral tract sampling. In the context 
of CT, a butterfly filter specific to CT machines was devel- 
oped to introduce a CT source[26]. Additionally, the program 
supports phase-space files for radiation field simulations and 
user-defined source parameters. Beyond these external ex- 
posure scenarios, the software offers two methodologies for 
internal exposure dose assessment: a swift computation ap- 
proach based on specific absorbed fraction (SAF) values de- 
termined via the Monte Carlo algorithm and a more detailed 
method reliant on the nuclide spectrum. 


3. Dose calculation algorithms 


The software offers three calculation modes: external, in- 
ternal, and environmental external radiation dose assessments 
[27]. Each mode facilitates computations of organ dose, ef- 
fective dose, and three-dimensional dose distribution. Organ 
doses can be directly calculated by utilizing voxels for Monte 
Carlo transport. Furthermore, to support three-dimensional 
radiation dose evaluations, THUDosepp offers both voxel 
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tally and mesh tally modes. In the voxel tally mode, when 
the input geometry is presented in voxels, the software gen- 
erates a three-dimensional matrix identical in size to the ge- 
ometry. It records physical quantities in each voxel grid and 
saves them in the respective element of the output matrix. 
This output matrix represents the three-dimensional statisti- 
cal outcomes of the required physical quantities. Addition- 
ally, when there is a discrepancy between geometric and tally 
resolutions, data can be documented in the mesh tally mode, 
which derives from the G4MeshKernelRect class. 

The calculations for organ and effective doses occur in the 
dose assessment module, termed as dose postprocessing. The 
aim here is to ascertain radiation-sensitive organ doses and 
then compute the effective dose, taking into account the tissue 
weight factor. Moreover, the software includes dose conver- 
sion coefficient computing capabilities, offering a theoretical 
foundation for radiation protection dose evaluations. 

Assessing the dose for the red bone marrow is both vi- 
tal and challenging in external radiation dose evaluations. 
THUDosepp adopts distinct methodologies for dose compu- 
tations for typical organs and the red bone marrow. For organs 
excluding the red bone marrow, an average dose methodology 
is employed, as detailed in equations (1-2): 


dé 

D = — 1 

7 (1) 

Dr= a / Ddm, (2) 
MT mr 


where dé denotes the mean energy imparted by ionizing 
radiation to matter of mass dm, D denotes the absorbed dose, 
and mr denotes the mass of an organ or tissue. These two 
equations were used to calculate the average absorbed dose of 
organ T, Dr. Furthermore, THUDosepp was used to directly 
count the energy deposition to obtain the absorbed dose of the 
organs with the exception of the red marrow and endosteum. 

For the red bone marrow, THUDosepp provides three 
evaluation methods. 


Direct method: The approach considers the average dose of 
bones containing red marrow as the representative dose 
for the red marrow itself, termed as the bone-specific 
weighted average dose. This is computed as per the 
following formula, equation (3): 


Erpm _ Epone X 'RBM 


(3) 


Dram = 


MRBM MRBM 


The variable Drgm represents the absorbed dose of 
the red bone marrow, while mrgm signifies the mass 
of the red bone marrow. Ergm and Epone denote the 
energy deposits in the red bone marrow and bone, re- 
spectively. Lastly, rRgm stands for the mass ratio of 
the red bone marrow to the bone. 
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3CFs-improved method: Given that the equilibrium condi- 
tion between the red bone marrow and trabecular bone 
is not met in phantoms, the S factor (or KS factor) is 
employed to calculate the dose of voxelized red bone 
marrow, termed the 3CFs method. Liu et al.[28] later 
refined this approach, resulting in the 3CFs-improved 
method. The absorbed dose Drmnp for the overall red 
bone marrow is the mass-weighted average of the red 
bone marrow dose at various locations. The associated 
formulas are presented in equation (4): 


where Di wp and D p4 denote the absorbed doses 
of the red bone marrow and spongiosa in bone i, the 
KS factor adopts the KS data of 44-year-old adults, Æ 
and ¢(£) denote the photon energy and fluence in the 
spongiosa, respectively. 


Drem = X RÌ - Dieu (5) 


1 


The final average dose to the total red marrow is pro- 
vided by equation (5), where R* denotes the ratio of 
the mass of red marrow in bone i to the total mass of 
red marrow. 


Dose-response function method: calculated with dose- 
response function method[29], as follows, equation 
(6-7): 


D(rr, x) = P(E, rs,x)R(rr + rs,z, E)dE, (6) 
E 


m(rr, x) 


m(rr) 


Dsralrr) = >> 


x 


D(rr, 2). (7) 


Bone-specific absorbed dose per particle fluence (Gy - 
cm?) to red bone marrow and endosteum is a function 
of particle energy and skeletal region, i.e., R(rr + 
rg,x,E) in equation (6). The fluence of particles 
of energy E in the source tissue rg of bone site x, 
®(E,rg,x), was simulated by THUDosepp. The ab- 
sorbed dose to tissue ry (red bone marrow) at bone 
site x was calculated using equation (6). Addition- 
ally, equation (7) is used to obtain the average absorbed 
dose of red marrow, where m(rr, x) and m(rr) denote 
the mass of tissue ry at bone site x and the total mass 
of rr, respectively. 
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For internal radiation dose computation, the software offers 
a quick estimation mode grounded on SAF values alongside 
a more detailed computation mode. The detailed approach 
utilizes Monte Carlo simulation, relying on the organ’s nu- 
clide spectrum. This strategy accommodates various models, 
working positions, and nuclide residence durations, employ- 
ing the organ nuclide spectrum for the Monte Carlo simu- 
lation. Subsequently, the software calculates the equivalent 
and effective doses committed to the organs irradiated within 
the phantom. However, given this aspect is beyond the main 
scope of this research, it will not be extensively discussed 
here. 


4. Three-dimensional visualization 


THUDose facilitates the visualization of geometry and par- 
ticle trajectories in the ’wrl’ format, compatible with most 
three-dimensional software. Moreover, TTHUDosepp en- 
hances the three-dimensional visualization module, offering 
windows in four distinct views: perspective, front, side, and 
top. This advanced view supports the customization of vis- 
ible slices and adjusts the displayed dose value ranges. For 
CAD geometries and polygon-mesh phantoms, 'THUDosepp 
allows imports in the ’ply’ format and offers the ability to 
display components in a variety of colors. When it comes 
to voxel phantoms, THUDosepp features a display func- 
tion for binary formats, where each organ can be distinc- 
tively showcased using a unique color set. In the realm of 
dose assessment, it is often pivotal to merge the geometry 
(or phantom) with the dose distribution, aiding in pinpoint- 
ing high-dose locations and facilitating radiation risk analy- 
ses. THUDosepp’s visualization module presents a solution 
that concurrently visualizes both geometry and dose distribu- 
tion. The integration of isodose lines further aids users in their 
analyses. Furthermore, while organ doses are often rendered 
as non-intuitive tables or datasets, THUDosepp revolution- 
izes this by offering three-dimensional visualizations of organ 
doses, enabling the display of specific organs to streamline 
user comparisons and in-depth analyses. 


C. Validation of THUDosepp 


The program’s accuracy in phantom dose assessment un- 
derwent rigorous testing in several scenarios. THUDosepp 
was benchmarked against various ICRP publications, cover- 
ing areas, such as external photon irradiation, external neu- 
tron irradiation, and internal electron irradiation, to authen- 
ticate its physical processes and dose calculation method- 
ologies. Additionally, by transforming the Chinese refer- 
ence polygon-mesh phantom into voxels, discrepancies be- 
tween the voxelized phantom and the established benchmark 
were scrutinized. Given the program’s reliance on a three- 
dimensional visualization platform, its prowess in depicting 
the geometry and results in 3D is also showcased. 
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1. Benchmark against ICRP publications 


The external irradiation benchmarking leveraged insights 
from ICRP publication 116. This evaluation contrasted dose 
conversion coefficients for ISO irradiation with monoener- 
getic photons and AP irradiation with monoenergetic neu- 
trons. The internal irradiation benchmarking drew from the 
research of Zankl et al.[30, 31], which examined the SAF 
values of organs when subjected to electron and photon ir- 
radiation. 

ICRP Publication 116 deduces reference conversion coef- 
ficients for both the effective dose and organ-absorbed doses 
under various external exposure conditions. This is based on 
the ICRP reference phantoms for both genders. It offers data 
from external beams of monoenergetic photons, electrons, 
positrons, and neutrons with configurations such as AP, PA, 
LLAT, RLAT, ROT, and ISO. In our assessment, the focus was 
on validating dose conversion coefficients specifically for the 
ICRP reference male, choosing monoenergetic photons and 
neutrons ranging from 10keV to 10GeV for the process. 

Zankl et al. employed the ICRP reference male and the 
EGSnrc Monte Carlo program to deduce SAF values for pho- 
tons and electrons, considering varied source organs and tar- 
get organs. To ascertain the precision of SAF values derived 
using THUDosepp, the ICRP reference male served as the 
basis for calculating the SAF values for electrons moving 
from the kidney to three distinct organs: the kidney itself, the 
adrenals, and the liver. The outcomes of this exercise were 
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2. Polygon-mesh phantom voxelization 


The polygon-mesh phantom was crafted by demarcating 
the organ boundaries of the voxel phantom using three- 
dimensional visualization software. This formed polygon 
mesh phantom boasts flexibility, allowing adjustments in 
height, weight, and posture. Once tailored, the modified 
polygon-mesh phantom was then imported into THU Dosepp 
for dose calculations by transforming it back into a voxel 
phantom. Illustrating this with the Chinese polygon mesh ref- 
erence male phantom (C RAMs) as a case in point, the steps 
for voxelization can be seen in Figure 3. 

First, the polygon-mesh phantom was adjusted using the 
three-dimensional modeling software Rhino to realize the tar- 
get size and posture. Layers were added and named for 
each organ. The Python plug-in of Rhino was used to write 
scripts for batch export based on the polygon mesh phantom 
in *.3dm form. Then, each organ was exported as an stl file 
with the same file name as the organ layer, and simultane- 
ously the stl file list was exported, i.e., x.lst text file, in which 
the file names are sorted from the largest to the smallest organ 
volume. 

Subsequently, preparation of a phantom information file 
is necessary. This file should detail the index and antic- 
ipated weight of each organ. Moreover, during the post- 
processing phase of dose calculations — particularly when 
merging residual tissues — it is essential to identify the group 
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Fig. 3. Flow chart of the voxelization operation in THU Dosepp. 


to which an organ belongs. This type of information is then 
saved in a text file, denoted as (* organmassOrg.txt). 

The previously prepared files were then loaded into 
THUDosepp. Users can set their preferred voxel size, ex- 
pressed in millimeters (mm). Once the conversion mode is 
chosen, the program permits a maximum permissible devia- 
tion from the intended organ weight. Itis worth noting that 
all organs are converted in descending order of size. Given 
this approach, conflicts may arise when smaller-volume or- 
gans intersect with larger ones. To address such conflicts, the 
program offers three resolution methods: overwriting, retain- 
ing, and randomizing. 


overwrite: ’Overwrite’ refers to the region where the con- 
flict occurs is covered by the small volume organ that 
is converted later. 


keep: ’Keep’ denotes that the current conversion result is 
maintained, and the small volume organ that will be 
converted later avoids the conflict area. 


random: ’Random’ means that in conflict areas, small- 
volume organs are randomly covered. 


The maximum deviation was used to ensure the upper limit 
of deviation between the voxelized organ mass and expected 
mass. The algorithm is realized by a loop that combines the 
corrosion and dilation algorithms. In this study, a polygon- 
mesh phantom C RAMs with a height of 170 cm and weight 
of 63 kg was voxelized. In the THUDosepp parameter con- 
figuration, the voxel resolution is 1.74lmm x 1.741mm x 
1.0007, the conversion mode is selected as ’overwrite,’ and 
the maximum deviation is set to 2%. The results of the output 
voxel phantom, CRAM, were compared for organ mass and 
tested for the efficiency of the entire conversion process. 


3. Three-dimensional dose assessment 


A three-dimensional depiction of an indoor work environ- 
ment was evaluated across several layers: the combined vi- 
sualization of ray traces with geometry, integrated display of 
phantom files with three-dimensional dose distribution, and 
independent portrayal of phantom organ doses. 
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For the ray trace and geometry visualization test, a new- 
born phantom undergoing radiography was used. This phan- 
tom was subjected to a vertebral photon beam in an AP con- 
figuration. The distance between the source and the phan- 
tom’s back was set at 100 cm, with a radiation field measur- 
ing 12.5cm x 17.4cm. To examine the integrated visualiza- 
tion of phantoms with three-dimensional dose distributions, 
and the distinct portrayal of organ doses, an external irradia- 
tion scenario was simulated using an isotropic point source. 
The chosen phantom for this test was a reference Chinese 
adult male. Positioned directly in front of the phantom was 
a '87Cs point source, emitting energy at 661keV. The three- 
dimensional dose distributions and organ doses were concur- 
rently scored and then utilized in the visualization tests. 


HI. RESULTS 


A. Benchmark against ICRP publication 116 


In the study, THUDosepp was employed to determine the 
absorbed dose of the ICRP reference male subjected to ISO 
geometry. We focused our analysis on radiation-sensitive or- 
gans with significant radiation weight factors. The dose cal- 
culation for the phantom was executed using 100 CPU cores, 
transporting 1E8 particles in a parallel manner. Figure 4 (a) 
displays the organ absorbed dose per fluence for 1E8 mo- 
noenergetic photons. Within this, both the ’Direct method’ 
and *3CFs-improved method’ were employed to assess the 
dose for the red bone marrow. These methods correspond to 
*THUDosepp’ and ’THUDosepp — 3CF’ in the legend, re- 
spectively. It is evident that our results align with the data 
from ICRP Publication 116. Notably, over 95% of our find- 
ings deviate from the ICRP Publication 116 data by less than 
5%. The remaining deviations, particularly in results with 
larger discrepancies, are predominantly attributed to substan- 
tial statistical errors in both datasets at lower energies. We 
further juxtaposed the results of red bone marrow dose calcu- 
lations using the 3CFs-improved method. This approach ef- 
fectively rectified the elevated values observed in direct cal- 
culation methods in lower-energy areas, as confirmed by a 
prior study[28]. 

Using analogous calculation parameters as applied for pho- 
tons, we determined the dose for neutrons irradiated with an 
AP geometry. Our focus was primarily on lung, skin, and 
brain tissues. Figure 4 (b) presents the absorbed dose per flu- 
ence for 1E8 monoenergetic neutrons. When juxtaposed with 
ICRP data, our findings exhibited consistency across a broad 
spectrum of energy levels, ranging from keV to GeV. 

ICRP Publication 133 offers data pertinent to internal ex- 
posure, including SAF values deduced via the Monte Carlo 
method. This represents the proportion of radiant energy orig- 
inating from the source region that is absorbed per mass in the 
target region. Give that these ICRP data stem from Zankl’s 
work, we reference them as ’Zankl et al.’ in our subsequent 
analyses. Figure 5 (a) and Figure 5 (b) illustrate the SAFs 
of electrons and photons traveling from the kidney to various 
target organs, including the kidney itself, the adrenal glands, 


202310.03382v1 


chinaXiv 


471 
472 
473 
474 
475 


476 


477 


478 
479 
480 
481 
482 


483 


(a) photon irradiation 


Energy (MeV) 


Energy (MeV) 


100 100} Breast 
10 10} 
1 4} 
0.1 0.1 
= — ICRP 116 vil. 
—— THUDosepp — ICRP 116 — ICRP 116 
om —— THUDosepp-3CF Sooty —— THUDosep, vet —+— THUDose,, 
1E-4 4. 1E-4 i. á - 1E-4 à 
0.01 0.1 1 10 100 1000 10000 0.01 0.1 10 100 1000 10000 0.01 0.1 1 10 100 1000 10000 


Energy (MeV) 


Ss 
= 


eS 
o 
= 


—— ICRP 116 


Sosti —— THUDose,, 


m 
E 


—— ICRP 116 
— +— THUDose,p 


| Remainder tissue 


— ICRP 116 
—— THUDose,p 


01 1 10 100 1000 10000 
Energy (MeV) 


0.01 


Absorbed dose per fluence (pGy-cm?) Absorbed dose per fluence (pGy-cm?) 


Absorbed dose per fluence (pGy-cm?) Absorbed dose per fluence (pGy-cm?) 


(b) neutron irradiation 


-æ 
o 
o 
© 


1000 


= 
o 
© 


100 


= 
G] 


10 


ICRP116 
= THUDOSepp 


Energy (MeV) 


10 100 1000 10000 0.1 1 10 100 1000 10000 


Energy (MeV) 


Absorbed dose per fluence (pGy-cm?) Absorbed dose per fluence (pGy-cm?) 


1000 


100 


10 


ICRP116 
= THUDose,, 


ICRP116 
= THUDose,, 


40 100 1000 10000 
Energy (MeV) 


Absorbed dose per fluence (pGy-cm2) 
Absorbed dose per fluence (pGy-cm2) 


Energy (MeV) 


0.01 0.1 1 10 


100 
Energy (MeV) 


1000 10000 


10 100 1000 10000 


Absorbed dose per fluence (pGy-cm2) 


Fig. 4. (a) Organ absorbed dose per fluence for mono-energetic photons incident in ISO, (b) Organ absorbed dose per fluence for mono- 


energetic neutron incident in AP. 


and the liver. When the kidney functioned as both the source 
and target, the relative deviation in SAF values was capped at 
3%. However, for target organs, such as the adrenal gland or 
liver, substantial deviations were observed at lower energies 
due to considerable statistical errors. Yet, the results closely 
matched at higher energy levels. 


B. Polygon-mesh phantom voxelization 


Upon voxelization of CRAMs, we obtained the voxel 
phantom CRAM with dimensions of [276, 156, 1700] and a 
cumulative count of 73195200 voxels. This encapsulates 103 
organs and tissues of the Chinese reference male. A com- 
parison of voxelization results with the ground truth values is 
depicted in Figure 6. The orange bars illustrate the genuine 
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value, which signifies the anticipated mass of the organ prior 
to voxelization. Conversely, the green bars indicate the post- 
conversion mass of the respective voxel, following the trans- 
formation from polygon-mesh organs. The relative deviation 
between these two, as showcased in Figure 6, aligns with the 
predetermined permissible deviation threshold of 2%. 

Furthermore, we chose four phantoms from the Chinese 
adult standing phantom repository for the voxelization pro- 
cess. Their respective conversion durations are detailed in 
Table 1. The process to voxelize an adult phantom, contain- 
ing between 50 to 90 million voxels, spans approximately 2-3 
minutes, contingent on the phantom’s size. 

The resultant visualization file is outputted in the *.wrl for- 
mat, allowing users to leverage any three-dimensional visual- 
ization software for analysis. This includes evaluating object 
coordinates, radiation field dimensions, and other pertinent 
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Fig. 5. (a) SAF values of electrons with kidney as the source organ, (b) SAF values of photons with kidney as the source organ. 


data. The geometry and ray trajectories of neonatal phan- 
tom radiographs focusing on the pelvis are presented in Fig- 
ure 7 (a). Within THUDosepp, the ’run parameter’ dictates 
whether the geometry is generated. If geometry generation is 
approved, then the default history of transported particles is 
set to ensure that the simulation scene can be rendered with 
minimal computational effort. The amalgamation and display 
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Fig. 6. Comparison of organ mass between voxelization results and 
truth values. 
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Table 1. Voxelizaiton cost time. 


Phantom Phantom size conversion 
time (s) 
H155M45 [258 146 1552] 120 
H165M67 [282 168 1650] 150 
H1I75M75 [288 184 1752] 156 
H185M8&4 [294 176 1850] 196 


of tracks and geometries significantly streamline project ex- 
aminations. 

Dose distribution within a phantom offers the most 
straightforward approach to three-dimensional dose evalu- 
ation. The amalgamation of the phantom and its three- 
dimensional dose representation facilitates this methodology. 
Figure 7 (b)-7 (c) illustrate the three-dimensional dose dis- 
tribution from a !37C's point source utilizing an AP geom- 
etry. Figure 7 (b) presents four distinct display modes of 
the phantom and dose: (i) showcases the phantom’s cross- 
sectional view, (ii) emphasizes the dose distribution, (iii) in- 
tegrates both phantom and dose visuals, and (iv) supplements 
(iii) with isodose lines. This versatile visualization technique 
assists users in analyzing radiation field metrics and dose lay- 
outs. A comprehensive three-dimensional perspective, en- 
compassing perspective, front, lateral, and top views comple- 
mented by a color scale, is portrayed in Figure 7 (c). 

Additionally, THUDosepp introduces a feature for visu- 
alizing the three-dimensional organ dose distribution. Figure 
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Fig. 7. (a) Geometry and track of neonatal pelvic X-ray photography, (b) Fused display of voxel phantom and three-dimensional dose 
distribution, (c) Three-dimensional dose distribution of isotropic point source irradiation with 137 C's, (d) Three-dimensional organ dose of a 


random energy spectrum. 


Š sos 7 (d) displays the organ doses across various energy spectra. 


s27 Each organ’s dose is represented by a distinct color, signify- 


= sæ ing its specific dose amount. 


IV. CONCLUSION 


529 


sso In this research, we introduced the Monte Carlo simula- 
531 tion platform, THUDosepp, tailored for three-dimensional 
s532 phantom dose evaluations. This platform is adept at parallel 
processing of photons, electrons, neutrons, and other parti- 
cle types using both CPU and GPU. It offers a plethora of 
methods for defining sources and geometries. Significantly, 
it incorporates a voxelization technique for the polygon mesh 
phantom, achieving phantom conversions with a precision un- 
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sxe der 2% in approximately 2 min. Moreover, it supports calcu- 
lations for both external and internal organ irradiation doses 
and enables visualization of three-dimensional geometry, par- 
ticle trajectories, and dose distributions. Catering to an array 
of requirements in three-dimensional phantom dose assess- 
ments, this platform proves invaluable for both radiation pro- 
tection and clinical scenarios. Currently, its applications span 
personal dose evaluations[32], environmental exposure dose 
coefficients[27], internal radiation cancer risk studies, dosi- 
metric reconstructions for a radiological mishap in Nanjing, 
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